Mechanisms of direct reactions with halo nuclei 
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Abstract. Halo nuclei are exotic nuclei which exhibit a strongly clusterised structure: they 
can be seen as one or two valence nucleons loosely bound to a core. Being observed at the ridge of 
the valley of stability, halo nuclei are studied mostly through reactions. In this contribution the 
reaction models most commonly used to analyse experimental data are reviewed and compared 
to one another. A reaction observable built on the ratio of two angular distributions is then 
presented. This ratio enables removing most of the sensitivity to the reaction mechanism, which 
emphasises the effects of nuclear structure on the reaction. 



1. Introduction 

The advent of radioactive-ion beams has enabled the exploration of the nuclear landscape far 
from stability. This technological breakthrough has led to the discovery of exotic nuclear 
structures, such as halo nuclei [1]. Halo nuclei exhibit a large matter radius in comparison 
with their isobars. This peculiarity is qualitatively explained by the small binding energy of one 
or two valence nucleons, which then have a significant probability of presence at a large distance 
from the core of the nucleus. They thus exhibit a strongly clusterised structure in which one or 
two loosely-bound nucleons form a diffuse halo around a dense core. Examples of one- neutron 
halo nuclei are ^^Be and ^^C, whereas ^He and ^^Li are the best known two-neutron halo nuclei. 
'^B is a candidate one-proton halo nucleus. 

Being far away from stability, halo nuclei cannot be studied through usual spectroscopic 
techniques, and we have to rely on indirect methods, such as reactions, to infer information about 
their structure. Breakup is probably the mostly used reaction to study halo nuclei [2, 3j. In that 
reaction, the halo dissociates from the core through the interaction with a target. In order to 
get valuable information from experimental data, various reaction models have been developed: 
the Continuum Discretised Coupled Channel model (CDCC) the Time-Dependent approach 
(TD) [5j, and the Dynamical Eikonal Approximation (DEA) [6]. These models are presented in 
Sec. [2j That section also includes a comparison between them for the breakup of ^^C on Pb at 
68AMeV [7], which has been measured at RIKEN [3]. 

To reduce the sensitivity of the measurements to the reaction mechanism, one can also look 
for an observable that is independent of the reaction mechanism. Such an observable is presented 
in Sec. [3l It corresponds to the ratio of two angular distributions. According to the prediction 



of the Recoil Excitation and Breakup model (REB) [8|, this ratio should depend only on the 
projectile structure. This prediction is confirmed by DEA calculations of reactions involving 
^^Be, the archetypical one-neutron halo nucleus [9]. 

2. Reaction models 

2.1. Theoretical framework 

Let us consider the reaction involving a one-neutron halo nucleus, in which a valence neutron / 
is loosely bound to a core c. Such a two-cluster projectile is described by the Hamiltonian 

Ho = Tr + Vcf{r), (1) 

where r is the neutron-core relative coordinate, Tr is the c-f kinetic-energy operator, and Vcj 
is a real potential that simulates the interaction between the halo neutron and the core. The 
parameters of this phenomenological potential are adjusted to reproduce the binding energy of 
the projectile and some of its low-lying levels. The structure of the target T is usually neglected 
and its interaction with the projectile components is simulated by optical potentials Vct and 
VfT- Within this framework, the theoretical study of reactions involving halo nuclei reduces to 
solving the three-body Schrodinger equation 

[Tr + Ho + Vct{Rct) + VfT{RfT)] ^{r,R) = Er^ir, R), (2) 

where R is the projectile-target relative coordinate and Tr is the corresponding kinetic-energy 
operator. Eq. ^ is solved with the initial condition that the projectile is in its ground state (po 

^(r,i^) ^ e^^^+-(/>o(r), (3) 

Z— ^— oo 

where K is the wave number of the initial relative motion of the projectile to the target, which 
is assumed to be in direction Z. The reaction models differ mostly in the way Eq. ([2]) is solved. 

In the CDCC model, Eq. ([2]) is solved by expanding upon the eigenstates of Hq [Ij. Since 
breakup leads to the dissociation of the projectile, a tractable description of the projectile 
continuum must be included in this expansion. That is why the continuum is discretised. This 
is usually achieved by dividing the continuum into small energy intervals called bins. Each bin 
wave function is obtained by averaging the wave functions of the continuum states within the 
corresponding energy interval. This leads to a set of coupled equations, which we solve using 
FRESCO [lOj. Apart from the continuum discretisation, the CDCC technique solves Eq. ([2]) 
exactly. It is a purely quantal model with no approximation on the projectile-target relative 
motion, and it does not imply any restriction on the beam energy. However, it is very expensive 
in a computational point of view, in particular at high beam energy. 

A simpler way to solve Eq. ^ is the Time-Dependent (TD) approach, which relies on a 
semiclassical approximation. In that approximation, the projectile-target relative motion is 
simulated by a classical trajectory R[t). Along that trajectory, the projectile is affected by the 
time- varying Coulomb and nuclear fields of the target. The evolution of its internal structure 
is then obtained by solving a time- dependent Schrodinger equation [5]. That equation must 
be solved for each possible trajectory, with the condition that the projectile is initially in 
its bound state. Various computer codes have been developed to numerically solve the time- 
dependent equation. In this study, we will use the one detailed in Ref. [5j. Thanks to its semi- 
classical approximation, the time-dependent method is faster than CDCC. However, due to this 
approximation, quantal interferences between the trajectories cannot be taken into account. 

Another way to avoid the heavy CDCC machinery is to use the eikonal approximation. At 
sufficiently high energy, one can assume that the projectile-target relative motion will be close 
to the initial plane wave ([3]). The main idea of the eikonal approximation is to factorise that 



plane wave out of the three-body wave function ^ and to assume that the remaining factor is 
smoothly varying with R. This leads to the Dynamical Eikonal Approximation (DEA), whose 
equation is mathematically equivalent to a time-dependent Schrodinger equation for straight- 
line trajectories [6j. As in the TD approach, the DEA equation must be solved for all the 
possible "trajectories". However, it remains a quantal approximation since it does not rely on 
any semiclassical approximation. This implies in particular that, contrarily to the TD approach, 
the DEA accounts for the quantal interferences between different "trajectories". However, since 
it relies on the eikonal approximation, the DEA is valid only at intermediate and high energies. 

2.2. ^^C Coulomb breakup 

To compare the three aforementioned models, we consider the breakup of the one-neutron halo 
nucleus ^^C on lead at 68AMeV [3]. The breakup cross section is plotted in Fig. [1] as a function 
of the relative energy E between the halo neutron and the ^^C core after dissociation. All three 
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Figure 1. Energy distribution for the breakup 
of on Pb at 68AMeV. Comparison of CDCC, 
TD, and DEA [7J. Data from Ref. [3] 



Figure 2. Angular distribution for the 
breakup of ^^C on Pb at 68AMeV. Compari- 
son of CDCC, TD, and DEA [7]. 



models are in excellent agreement with each other [7] and with experiment [3]. This result shows 
that the reaction process is well understood, and that ^^C exhibits a two-cluster structure. 

Fig. [2] shows the breakup angular distribution, i.e. the breakup cross section as a function of 
the scattering angle 6 of the ^^C-n centre of mass [7]. CDCC and DEA lead to very similar cross 
sections. In particular they exhibit the same interference pattern, which confirms the validity 
of the DEA at these intermediate energies, even though it is much less time consuming than 
CDCC. As expected, the TD technique does not exhibit any quantal interference. Nevertheless, 
it reproduces the general trend of the quantal angular distributions, which explains why, once 
integrated over 9, it gives reliable results. 



3. Ratio technique 

Albeit reliable, the aforementioned models remain sensitive to their input parameters. In 
particular, the choice the optical potentials Vct and Vjt can affect significantly the results 
of the calculations [11] . An observable less sensitive to the reaction mechanism would naturally 
eliminate such sensitivity and improve the analysis of exotic nuclear structures from reactions. 

It has been observed that angular distributions for elastic scattering and breakup exhibit very 
similar features: Coulomb rainbow, oscillatory patterns etc. [12{ . This result can be qualitatively 
understood within the Recoil Excitation and Breakup model (REB) [8]. That model relies 
on two simplifying assumptions. First it includes an adiabatic treatment of the excitation of 
the projectile, i.e. the projectile excitation energy is neglected compared to the beam energy. 



Second, it neglects the interaction between the halo neutron and the target Vjt- Under these 
two assumptions an elegant factorisation for the elastic-scattering cross section is obtained [8] 



which is the product of a cross section for the elastic scattering for a pointlike projectile 
{da/dQ)p^^ times a form factor, which accounts for the halo extension 

i^o,o = / \Mr)\'e^^-''dr, (5) 

where Q is proportional to the transferred momentum. A similar factorisation can be performed 
for the angular distribution for the breakup of the projectile to energy E in the continuum 

where the form factor reads 



dEdn ' \dnj 



l^rf , (6) 
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^ijUE,r)Mr)e'^-''dr 



(7) 



It contains, besides the wave function of the initial bound state of the projectile (po, the wave 
functions describing the broken up projectile in partial wave Ijm at energy E in the continuum. 

The factorisations dH) and Q explain the results of Ref. [12] since both expressions contain 
the same {da/dO,)^^, which is responsible for most of their angular dependence. In addition, the 
REB model also provides the idea for the ratio technique: within that model, the ratio of cross 
sections ^ and ^ is the ratio of form factors ([5]) and ([7]), which depend only on the projectile 
structure. Accordingly, the ratio should be independent of the reaction process. In order to 
test this idea, we perform reaction calculations within the DEA [6], which does not include 
an adiabatic treatment of the projectile excitation, and which includes VfT- As a testcase, 
we consider the collision of ^^Be with Pb at 69AMeV, which correspond to the conditions of 
a RIKEN experiment [2] with which the DEA calculations are in excellent agreement . 
The result of this test is summarised in Fig. [21 which presents the DEA angular distributions, 
their ratio, and its REB prediction [U]. Note that instead of the mere elastic-scattering cross 
section ([3D, we use the summed cross section dasum/di^, which corresponds to all quasi-elastic 
processes (elastic and inelastic scattering and elastic breakup) . It is plotted here as a ratio to the 
Rutherford cross section (dotted line). This gives the best agreement between DEA and REB [9J. 
As already noted in Ref. [T2|, both angular distributions exhibit similar features. Interestingly 
their ratio (solid line) does not exhibit such features and it behaves smoothly with the scattering 
angle 6. In addition the DEA ratio is in excellent agreement with its REB prediction (thick 
grey line). This result suggests the ratio idea to be an excellent reaction observable to study 
loosely-bound systems, such as halo nuclei. Further tests have confirmed its independence to 
the reaction mechanism: the ratio computed for a carbon target is nearly superimposed to the 
one obtained on lead |9J. Moreover, the ratio is strongly sensitive to the binding energy, orbital, 
and radial wave function of the valence neutron [9], indicating that valuable information about 
the projectile structure can be inferred from this observable. 



4. Summary and prospect 

Halo nuclei are among the most peculiar nuclear structures discovered thanks to the development 
of radioactive-ion beams. They exhibit a strongly clusterised structure, which explains their 
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Figure 3. DEA calculations of ^^Be impinging on Pb at 69AMeV. Angular distribution for 
breakup at ii^ = 0.1 MeV (dashed line) and the summed cross section (dotted line) are shown 
with their ratio (solid line) and the REB prediction (thick gray line) [9]. 

large matter radius in comparison to their isobars. Being located on the shoulder of the valley 
of stability, they are studied mostly through reactions, such as elastic scattering and breakup. 
Therefore, accurate reaction models must be available to extract valuable information from 
reaction data. In this contribution the mostly used reaction models have been presented and 
compared for the breakup of ^^C, a typical halo nucleus, at intermediate energy: the Continuum 
Discretised Coupled Channel model (CDCC) [1], the Time-Dependent approach (TD) [5], and 
the Dynamical Eikonal Approximation (DEA) f6]. Although based on different assumptions, 
they lead to nearly identical cross sections, but for the quantal interferences, which cannot be 
reproduced within TD, due to its relying on the semiclassical approximation. 

To remove the influence of the reaction mechanism on observables, we suggest to measure the 
ratio of angular distributions, e.g. for scattering and breakup. As predicted by the REB model 
[8], this new observable is nearly independent of the reaction mechanism and provides detailed 
information on the halo structure. It seems that the idea could be easily applied to Borromean 
neutral systems like ^^Li. Extensions to charged valence nucleon, such as proton haloes may also 
be possible. We believe this ratio method will open a new era in the study of exotic systems. 
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